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Abstract 

We evaluate the fermion propagator in parity-conserving QED^ with 
N flavours, in the context of an IR domain approximation. This pro- 
vides results which are non-perturbative in the loopwise expansion 
sense. We include fermion-loop effects, and show that they are rele- 
vant to the chiral symmetry breaking phenomenon, that can be un- 
derstood in this context. 

Parity-conserving QED in 2 + 1 dimensions, in the context of the large- 
approximation, is an effective theory description for many interesting models 
of condensed matter theory p. Moreover, it provides a simple yet non-trivial 
model where to study chiral symmetry breaking. Unlike in 1 + 1 dimensions, 
there is a truly dynamical gauge field; besides, the model shares some im- 
portant properties with its 3 + 1 dimensional counterpart. For example, the 
dimensionful coupling constant for the 2 + 1 dimensional model plays an 
analogous role to the renormalization scale of the 3 + 1 dimensional theory. 
It is, however, a much simpler framework where to study chiral symmetry 
breaking, and it allows for the use of controlled approximations. Many de- 
tailed studies have been carried out using the truncated Schwinger-Dyson 
equations for the gap function S, both at zero [21- 11] and finite 0-^01 tem- 
perature. Results (which are consistent with the continuum calculations) 
have also been obtained by using numerical lattice calculations pT]-|15j. 

*0n leave from Centro Atomico Barilochc, 8400 S. C. dc Barilochc, Argentina. 
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It is well known that QED^ has a non trivial IR structure, where the 
a priori severe IR divergences can be consistently treated thanks to the 
model superrenormalizability ^Hj. Indeed, the leading and first sub-leading 
logarithms have been summed to all orders for the fermion self-energy |17j . 

The relevant degrees of freedom driving the chiral phase transition are the 
low-momentum ones, with an effective dynamical cutoff set up by the (square 
of the) coupling constant. It then makes sense to attempt approximations 
based on that property. In this note, we calculate the full fermion propa- 
gator in parity-conserving large- QED^, in the context of the IR domain 
technique presented in jTH], and comment on the relevance of this result for 
the chiral phase transition. 

The model may be defined in terms of an Euclidean generating functional 
Z{j,f],r]), which in the path-integral representation is given by 



Z{j,r],ri) = J [DA]V^Vtjj exp 



where the Euclidean action 5" is, 

S{i>,^lj,A) = + SgiA) (2) 

with Sf and Sg denoting the fermion and gauge field actions, respectively. 
For Sf, we assume a Dirac-like form: 

N 

Sf{i/j,^lj,A) = d'xJ2M^ + 'er,4)i^a (3) 

^ a=l 

while the gauge field has a standard Maxwell action: 



Sg{A) = d'x- F^^F^, , F^, = d,A, - d,A^ . (4) 



1 
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Here ejv = eA^~^/^, where e is a coupling constant (with the dimensions of the 
square root of a mass). It is the only dimensionful parameter in the model. 
The N~^/'^ factor is introduced, as usual in the large- expansion, in order 
to have a finite limit for the leading term in the expansion, when iV — >■ oo. 

The Dirac 7-matrices are in a (reducible) 4x4 representation of the Dirac 
algebra: 

= ( ? 4 ) ' = ( -a, ) ' = ( ? -a, ) 
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where cTj, j = 1, 2, 3 are the standard Pauh matrices. In this representation, 
the model has an extra 'chiral' symmetry, generated by a '75' matrix 

f 12X2 \ 
= [ 12X2 j ' 

which together with other generators complete the Lie algebra of an inter- 
nal U{2N) symmetry group. This symmetry is known to be dynamically 
broken|2j, and a non- vanishing mass for the fermions is generated as a con- 
sequence. The U(2N) symmetry breaks down to U{N) x U{N). 

Besides, a fermion mass does not break the parity symmetry of the 
massless model. This means that parity symmetry is not anomalous, since 
there are regularizations that indeed preserve that symmetry, like the Pauli- 
Villars method, for example. This guarantees that there is no radiatively 
induced Chern-Simons term, and that it is perfectly consistent to use a 
parity-conserving action, since parity is not broken by renormalization ef- 
fects. 

The symbol [2^^/^] in O denotes the gauge field functional integration 
measure, plus any gauge-fixing artifacts. An effect of the breaking of chiral 
symmetry is given by the generation of a non-vanishing mass gap S for the 
fermionic field. This may in turn be observed by evaluating the full fermion 
propagator G in the theory defined by ((H) 

G{x -y) = {^l^{x)i,{y)) , (7) 

which, (by integrating out the fermion fields) may be written as 

G{x-y) = [VA] {xW-'\y) (det exp[-5,(A)] , (8) 

where 

M = j [VA] (det eM-SM)] , (9) 

and we have adopted Schwinger's notation {x\K\y) for the coordinate space 
matrix elements (or kernel) of an operator K\ in the case above the fermion 
propagator in an external gauge field. 

In the large- iV limit, the leading contribution from the fermionic deter- 
minant comes from its quadratic part, namely, 

(det^)^ = exp(iVTrln^) = e^^[-W^'^\A) + 0{N-^)] (10) 

where 

W^'\A) = -^Tr[r' 4r'4\ . (11) 
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The exact form of the quadratic functional W^'^\A) is well known; it may be 
written as follows: 



W^'\A) = d'xF,^-^F,^ (12) 



The operator ^ , which has a nonlocal kernel in coordinate space, may be 
thought of as being defined by its Fourier representation: 

where it is obviously local. 

Thus, in the large- iV limit, the gauge field average of (jH)) has a purely 
Gaussian weight. The evaluation of the gauge field average may be further 
simplified by using a convenient gauge-fixing term: 

[DA] = VA exp[-l jd'xd-A{l + ^-^) d ■ A] , (14) 

where ^ is a real constant. Then 

G{x-y) = VA{x\^-'\y) exp[-SG{A)] , (15) 

where the quadratic functional Sg^A) may be written explicitly as 

Sg{A) = d?xA^M^,A, (16) 



with 



and 



M^, = M(v^) 



(17) 



e2 



M(V^) = (V^ + - 



In order to evaluate G, in particular its self-energy part, we follow the 
approach of jTHj, although with a difference that turns out to be crucial 
in our case, namely, keeping the fermion determinant contribution to the 
effective gauge field action. We remark that it is totally justified to ignore the 
determinant in ^Hj , since the fermions are assumed to have a non- vanishing 
renormalized mass, and that mass sets up a scale where to make the relevant 
IR approximations. In our case, however, the only scale appears in the one- 
loop contribution. Moreover, there will also be a non- vanishing fermion mass. 
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but here it will appear in a self-consistent way. This self-consistency depends 
upon the inclusion of the full one-loop gauge field propagator. Indeed, this is 
known to be important in any study of chiral symmetry breaking in QED3 j2] , 
and this calculation is no exception, as it will become clear in what follows. 

Following ^Hj, we represent the fermionic propagator in an external field 
in terms of a path integral over trajectories, as follows 



(20) 



I dT I VxVp exp[z / drp ■ x] V exp[—i / (ir^(r)] 
Jo+ J Jo Jo 

X exp[-ieN / drx ■ A{x{t))] , (19) 
Jo 

where V denotes path ordering, and the functional integral over is taken 
over paths with the boundary condition x^(0) = and x^{T) = y^. There 
is no boundary condition for the values of p^{t) at and T. It is clear that 
the functional average over A will only affect the last factor in p9j) . Indeed, 

/•oo f i*T 

G{x — y) = dT VxVp exp[i / dTp{T)-x{T)] 
Jo+ J Jo 

r 

X V exp[—i / dr ^(r)] exp(— JF) , 
Jo 

where ^ is a functional of the trajectories x{t) and a function of T, and is 
defined by the average: 

exp{— JF} = (^exp[—ieN J drx ■ A{x{t))]^ , (21) 

with 

^ JVA...exp[-SaiA)] 
^•••^ JVAexp[-Sa{A)] ' 

Being Sq quadratic, the integral over A can of course be performed exactly. 
The result for JF may then be written as follows 

= ^ljd\j d\,j,{u){u\M;,'\v)Uv) (23) 

where we have introduced the current associated to the particle's worldline: 

j^{u) = I dTX^{T)5[u-x{T)], (24) 
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with = Taking into account the fact that the worldhne begins at y 
and ends at x, we see that 

d^JM = 6{u- y) - 6{u- x) , (25) 

and because of this, we may evaluate the contribution of due to the longi- 
tudinal part of {M~^)fj_i, exactly. Namely, the longitudinal part of the propa- 
gator, which is where all the gauge-dependence is contained, will contribute 
to a term, which we denote J-'i, in the following expression for JF: 



where 



^ = + (26) 
jr^ = ej^ I dr [ dr' x{t) ■ x{t') {x{T)\A'r^\x{r')) (27) 







and 

= el (e 



:28) 



Each separate term in ()28p. when written in momentum space, is divergent. 
However, their sum may be evaluated by combining the two integrands under 
a common momentum integral, what yields a convergent result. The reason 
is, of course, that the first term in exactly kills the logarithmic divergent 
of the second. 

The momentum integration can then be performed, obtaining: 

m = ;^(e-l)//(r), (29) 

where 

m = (7-l + ln(^) + ^[l-cos(^)] 

+ " sin(— )ci(— )] I , (30) 

where 7 is the Euler's constant, r = |x — ?/|, and si (ci) denotes the integral 
sine (cosine) function. 

Since Ti does not depend on the trajectories x{t\ it is clear that the 
only contribution relevant to the fermion mass will come from Tt in fl27|) . To 
calculate JF^, we follow the procedure described in ^H]. After the approxima- 
tions are implemented, and a 'ribbon' regularization with a small-distance 
cutoff h is introduced [IHj, we end up with 

Tt = ^ £ dr dr' M"\y/\r - rf + 6^) (31) 
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where M ^{x) = {x\M ^|0). The kernel M(V— <9^) in coordinate space can 
be computed, yielding 



2 2 2 2 2 ■ 

2Cos(— ) + ci(— )sin(— ) - si(— )cos(— ) 



(32) 



where r = It is evident that, if we use the exact expression for M~^{x) 
in (jHT|) , we cannot evaluate the integral over r and r' analytically, something 
that we do need in order to obtain the fermion self-energy. One can easily 
check that the small-r behaviour of (j32|l corresponds to the case of a pure 
Maxwell gauge field action, namely, 

M-\x) ~ ^, r^O, (33) 

while it yields the pure one-loop action dependence for large-r: 

M-\x) r»e-2. (34) 

It is possible to use a simple analytic approximation to (jS21), which however 
captures both ()33|) and ()34|) exactly and does not differ significantly from 
the exact function in between. We use an approximant to M~^(x), which we 
denote by Ai~^{x), and is given by the expression: 

M^\x) = 1 . (35) 

Inserting this expression in place of the exact one, in (jS^, one can 

evaluate exactly the integral over r and r'. Of course, not all the dependence 
in h is interesting, since, in the end, we will renormalize and take the limit 
6 — >■ 0. Thus we may just collect the leading terms when T » b. Then 



Attn 



.32 1 8 ne^T, 
2^'^y~[^> 



(36) 



where terms that tend to zero when b/T —>■ have been neglected. It is then 
possible to insert (jHBj) and (j2HI) into (pUj) . and integrate over T, to obtain: 

^ r 4 
G{x-y) = k{j)^n exp^---^{^-l)fi{r] 

"^'^ e'^<^-y^ (z ^+S)-(i+^) , (37) 



(2vr) 
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where 

S^jf-ta^i-). (38) 

and K is a numerical constant. 

It should be evident from the above that S, modulo a renormalization 
that we shall explain in what follows, signals the existence of a mass for 
the fermionic field. It may be argued that a fermion mass appears just 
because the ribbon regularization breaks chiral symmetry. However, this 
only explains the existence of a perturbative mass, and in fact can (and 
will) be taken into account by a careful choice of the finite counterterms. 
The situation is analogous to the one in QED^ with a non gauge-invariant 
regulator: in general, a mass counterterm is usually required in order to have 
a gauge invariant renormalized theory. 

On the other hand, a non-zero (renormalized) mass which is non-perturbative 
in the coupling constant shall be an a posteriori justification for the use of 
this IR approach. 

S determines the position of the singularity in momentum space, and it 
is evident from (jHHj) that it needs a renormalization, since it diverges when 
6 — s> 0. The renormalized mass, which we shall denote by m, will differ from 
S by a subtraction corresponding to a mass counterterm: 

m = S-So (39) 

where Eg is the counterterm. Now such a divergence in this superrenormal- 
izable theory does only occur for the self-energy diagram which is of order 
e'^/N. Noting on the other hand that S may be rewritten as 

, 1 , e2 ,32 /i 

ln( — ) H ln( — 

47iN ^bfi' 47iN ^ 71 



^ = TZTr H—) + -TZTr ln(^^) (40) 



where /i is a mass scale, we see that the divergent part ^^ln(^) is indeed 

2 

of order ^, as expected in this superrenormalizable theory. Then we fix the 
divergent part of the mass counterterm to be such that 



m = niQ + -j-^ ln(— -^) . (42) 



so that 

32 /i 

where mo is a finite counterterm, which has to be fixed by a normalization 
condition. It should of course be of the form: 

'"o—C-. (43) 
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o .ij^ lib -r uiib / A A\ 

m = mo + - — -7 m( ) . (44) 



where C is a number of order 1. We do not want to introduce an extra mass, 
scale into the theory, so we have to use fi = m + 5m as the renormahzation 
scale: 

32 m + 5m, 

The reason why /x cannot be just equal to m is the following: m = has 
to be one of the possible solutions to equation (jH)), so that chiral symmetry 
is not explicitly broken. But setting m = in certainly requires a non 
vanishing 5m: 

5m = ^ exp(-47riV^) , (45) 

what determines the relation between 5m and mg. This is in fact the renor- 
mahzation condition that guarantees that the chiral symmetry is not explic- 
itly broken by the renormahzation procedure. Using relation we may 
finally write the expression for m as a transcendental equation: 



m = e 

32 



Tfi 

exp(47rAr— ) - 1 . (46) 



e2- 



where C cannot be determined within the approximation made in obtaining 
()3ip . Indeed, in this approximation the induced mass is a constant rather 
than a momentum dependent function. Therefore it impossible to decide 
which is the value of the constant C, related to the dynamical cutoff (~ |^), 
which is observed in the Schwinger-Dyson equations approach jH]. 

However, from a simple analysis of we see that besides the trivial 
solution m = 0, there will be a non-vanishing mass as long as remains 
smaller than a critical value iVc, given by: 

iVc = ^e^. (47) 

This result is consistent with the ones obtained in more standard settings; 
for example, the value Nc = fl of [Sj, corresponds to C = ln4 ~ 1.39. It 
is interesting to compare this result with the one one would have obtained 
by ignoring the one-loop correction to the photon propagator. As already 
mentioned above, there is no induced mass term in such an approximation, 
since one finds no term linear in T in ()36|1 . 

We conclude by stressing that this IR approximation scheme looks promis- 
ing as a tool to study the chiral phase transition in QED^ like models, in 
a gauge invariant setting. However, in order to obtain the full momentum 
dependence of the fermion mass, a higher order approximation than the one 
used here may be required. Indeed, it is plausible that fermion recoil effects 
are necessary to sense that momentum dependence. This effect can not be 
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appreciated in the approximation used here since the motion of the fermions 
is not perturbed by their interaction with photons [TH] . 
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